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We study a codimension 2 braneworld in the Einstein Gauss-Bonnet gravity. We carefully examine 
the structure of possible singularities in the system which characterize the braneworld through 
matching conditions. Consequently, we find that the thickness of the brane can be incorporated as 
the distributional source, which we dub quasi-thickness. On the basis of our formalism, we analyze 
the linearized gravity and show the conventional Einstein gravity can be recovered on the brane. In 
the nonlinear regime, however, we find corrections due to the thickness and the bulk geometry. We 
also point out a possibility that the thickness plays a role of the dark energy/dark matter in the 
universe. 
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I. INTRODUCTION 

The old idea that our universe may be a braneworld 
embedded in a higher dimensional spacetime is re- 
newed by the recent development in string theory 
which can be consistently formulated only in 10 dimen- 
sions Instead of 10 dimensions, however, most stud- 
ies of braneworld cosmology have been devoted to 5- 
dimensional models |3|]. This is due to the difficulty in 
treating the higher codimension object in a relativistic 
manner 4]. Nevertheless, it is important to explore the 
possibility of the higher codimension braneworld [5(- 

It is the Einstein Gauss-Bonnet gravity that is a nat- 
ural framework in 6-dimensions. Because a natural ex- 
tension of Einstein gravity to higher dimensions is the 
Lovelock gravity and it reduces to the Einstein Gauss- 
Bonnet gravity in 6-dimensions [(J. Moreover, the Gauss- 
Bonnet term appears in the low energy effective action 
in string theory. Therefore, it seems reasonable to con- 
struct a codimension 2 braneworld in the Einstein Gauss- 
Bonnet gravity. In fact, Bostock et al. claimed that 
6-dimensional Einstein Gauss-Bonnet gravity leads to a 
thin braneworld where a conventional Einstein gravity 
holds . More interestingly, they suggested that the de- 
viation to the conventional Einstein theory can be possi- 
ble when the thickness of the brane is taken into account 
and consequently the variation of the deficit angle is al- 
lowed. However, they have never given a precise scheme 
to calculate these corrections. Rather, they commented 
that the corrections can not be fixed exactly by the bulk 
equations. Therefore, it is important to construct a vi- 
able concrete model for the thick braneworld and clarify 
if corrections can be obtainable or not. 

In this paper, we investigate a codimension 2 
braneworld in the Einstein Gauss-Bonnet gravity tak- 
ing into account the thickness. It is difficult to treat 
the finite thickness as it is. However, by examining the 
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structure of the singularity in the equations of motion, 
we find a possibility to treat the thickness within the 
context of the distributional source. We name it a quasi- 
thick braneworld. This is the important point in our work 
which gives a framework to clarify various issues in the 
codimension 2 braneworld. In the case of the linearized 
gravity, we clarify the effect of the bulk on the brane by 
solving the whole set of equations of motion in the bulk 
with proper considerations of the boundary conditions. 
It turns out that the conventional Einstein gravity is re- 
covered at the linear level. However, we also show some 
corrections due to the thickness can be expected at the 
second order level. The effect of the bulk geometry in the 
nonlinear regime is also discussed. We conclude that cor- 
rections to the conventional Einstein theory can be fully 
obtainable in the case of the quasi-thick braneworld. We 
also point out an interesting possibility that the thick- 
ness plays a role of the dark energy/dark matter in the 
universe. 

The organization of this paper is as follows. In sec. II, 
we present a model for thick braneworld and derive nec- 
essary equations. The structure of the singularity is 
carefully examined. Then, we introduce an idea of the 
thick braneworld within the context of the distributional 
source, which we dub quasi-thickness. In sec. Ill, we 
set the background vacuum spacetime. In sec. IV, we 
solve the bulk geometry using linear perturbation the- 
ory. In sec.V, the effective theory for the codimension 2 
braneworld is presented both in the linear and nonlinear 
regime. The consistency of our formalism is emphasized 
here. In the final section, we summarize our results and 
discuss possible applications and extensions. 



II. QUASI-THICK BRANEWORLD: A MODEL 
FOR THICK BRANEWORLD 



We consider a codimension 2 braneworld with a pos- 
itive tension in the 6-dimcnsional bulk spacetime. The 
Lovelock gravity is a natural framework in higher dimen- 
sions and reduces to the Einstein Gauss-Bonnet gravity 
in 6-dimensions. Moreover, the Gauss-Bonnet term ubiq- 
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uitously appears in the low energy limit of string theory. 
For this reason, the Gauss-Bonnet (GB) term is intro- 
duced in our model which is described by the action 



-g a 



(i) 



where k 2 is the 6-dimensional gravitational constant, 
3(6W and g^ are the 6-dimensional bulk and the our 
4-dimensional brane metrics, respectively. Here, £ ma tter 
is the Lagrangian density of the matter on the brane, and 
a is the brane tension. The GB term is given by 



r>2 rya 

H GB — H 



bed 



Rabcd — 4i?° R a b 



R 2 



(2) 



The Latin indices {a, &,•••} and the Greek indices 
{/i, ■ ■ ■ } are used for tensors defined in the bulk and 
on the brane, respectively. The 6-dimensional Einstein 
equation derived by varying the above action with re- 
spect to g^ takes the form 

G a b + aH ab = n 2 T ab , (3) 

where a is the GB coupling constant with dimension 
[a] = L 2 and 



:9abRc!B + %RRab — 4i? a rfi? d 6 



-4R ae R a dbe + 2R a d ' ° Rbdef 



(4) 



is an analogue of the Einstein tensor stemmed from the 
GB term ©. 

We will assume that a 6-dimensional metric has axial 
symmetry, which reads 



ds 2 = dr 2 + g» v {r, x^dx^dx" + L 2 (r, x^dO 2 



(5) 



This assumption corresponds to the Z2 symmetry in the 
Randall-Sundrum braneworld model. Here we have in- 
troduced polar coordinates (r, 9) for the two extra spatial 
dimensions, where < r < 00 and < 9 < 2ir. As we lo- 
cate a four-dimensional brane at r = 0, which is a string 
like defect, we must take the boundary condition 



limX(r,i")=0, lim L'(r, x») 



const. 



(0) 



where the prime denotes derivatives with respect to r. 
The first condition realizes the 4-dimensional brane at 
r = 0. And the second condition allows the existence of 
the conical singularity. The structure of conical singular- 
ity is a 2-dimensional delta function 5(r)/L. As we will 
sec later, the deficit angle is determined by the tension 
of the brane. 

Possible components of the energy-momentum tensor 
T a0l which could be balanced with the singular part of 
Einstein tensor, are given by 



T ab (singular) = ( T„ v (singular) ) .(7) 

T ee (singular) 



Notice that (r, r) component does not appear in the 
energy-momentum tensor Q, because this component 
must be balanced with the (r, r) part of Einstein tensor 
which consists of only the first derivatives with respect 
to r and hence there is no chance to have a singularity. 
On the other hand, (9, 9) component may have a singu- 
larity. To determine the structure of the singularity in 
the energy-momentum tensor, we need to examine the 
structure of the singularity in the Einstein tensor. The 
singular part of the Einstein tensor gives the matching 
conditions. Using the metric the non- linear matching 
condition for [9, 9) component is written as 



-K' + a 



^K a pK^IC a - 2KK l3 a K a p + ^K 3 



ARP a R a p - 2KR 



k 2 T\ (singular) ,(8) 



where K is the trace of is the extrinsic curvature, K^ v = 
— l/2g^ v ^ r . The (//, v) components of matching condi- 
tions read 
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K» v - 5»K 



-4a 



„ DL Z>i 



-2a K^^pK^a - 2K^ a K a f jK u + 2KK^ a K a v 

-K 2 K^ V - 2K a fj R^ va + 2KR fl v 
-2K^ a R a v - 2K a v R^ a + K» V R 



2K^ a R a p - KR 



K 2 T^ (singular) , (9) 



where | denotes the 4-dimensional covariant derivative 
and W^v is defined as 

W» v = K» a K a u - KK» V + (K 2 - K a p K p a ) (10) 

In general, A" M „ and L'/L could have discontinuities at 
r = 0, namely, limA^r — e,x^) ^ K^r — 0,x^) and 

e— >0 

the same for L'/L. For this case, two kinds of singular 
structures exist in Eqs. JSJ) and @ under the boundary 
condition ©■ One is a 2-dimensional delta function sin- 
gularity 8{r)/2-KL like K'/L, L"/L, etc. The other is a 
1-dimensional singularity S(r) such as K', which is less 
singular. These kinds of singularities require the same 
kinds of singular contribution in the energy-momentum 
tensor. The origin of two kinds of singularities in the 
energy-momentum tensor can be understood by the fol- 
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lowing schematic Taylor expansion: 



= e,x^ 



1 (o) 

-T ab (r = 0,1") 



K 1 ) 

-T ab {r = 0,x») , (11) 



where T ab is the energy-momentum tensor of smooth 

(o) (i) 

matter, and T ab and T ab denote the Taylor coefficients. 
Here, e represents the thickness of the matter distribu- 
tion. In the conventional thin limit, only the first term 
is usually taken. We propose to incorporate the second 
term to take into account the thickness of braneworld. 
Thus the possible energy-momentum tensor takes the 
form 







Tab (singular) 



T, 




where represents the conventional brane matter, 
while S^v and Sgg describe the extra matter which mim- 
ics the thickness of the braneworld. Note that Tgg can 
not be allowed due to the structure of Eq. (jHJ) . It should 
be stressed that this ansatz completely makes sense in 
the axially symmetric spacetime. In this way, the thick- 
ness can be treated within the context of a distributional 
source, which we dub quasi-thickness. In a sense, we are 
considering the "internal structure" of the braneworld. 



III. VACUUM BRANEWORLD 

The linear analysis is a first step to understand the ef- 
fective theory on the braneworld. As a background space- 
time, let us first consider the vacuum braneworld. Later, 
we will consider the fluctuations around this vacuum so- 
lution. Then, we can proceed to the discussion of the 
nonlinear dynamics of the codimension 2 braneworld. 

Because of the symmetry of the vacuum, the metric 
can be expressed as 



ds 2 



dr 2 + a 2 {r)r] llv dx' 1 dx v + b 2 {r)d6 2 , (13) 



where a(r) and b(r) depend only on r. The energy- 
momentum tensor (|12|) of the vacuum brane can be char- 
acterized by the tension a: 



T a b (singular) 





v 2tto 





(14) 



Here, we have introduced only the 2-dimcnsional distri- 
bution. This seems a reasonable assumption because the 
vacuum can not have an internal structure. 

Given the metric (|13fl . the Einstein Gauss-Bonnet 



equations ©, off the brane, in the bulk are written as 
2 
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a 2 b 

a' 3 b' 
~aW 



. a" a' 2 







The solution is obtained as 

ds 2 = dr 2 + rj^dx^dx" + c 2 r 2 d6 2 



,(15) 
, (16) 
0. (17) 

(18) 



where c is a constant of integration. For c ^ 1, we have 
a conical singularity at r — 0. This can be seen in the 
following way: First, we specify the boundary condition 
at r = by following the standard procedure, b'(r = 0) = 
1. Then, the 2-dimensional curvature is calculated as 



b 



,. 1 fb'(r = e)-b'(r = 0) 
lim — 1 

e^O b 



(c-1) 



6(r) 



(19) 



Now, the matching condition (|5J) determines the deficit 
angle in terms of the brane tension as 



1 



2tt 



(20) 



Note that c < 1, because the tension a is positive. This 
result gives the insight that the deficit might depend on 
when the inhomogeneous matter exists on the brane. 



IV. LINEAR PERTURBATION ANALYSIS 

Having considered a background solution and found 
that the deficit angle seems to depend on the brane mat- 
ter, we now consider a perturbation around this solution 
to examine the behavior of the deficit angle fluctuations. 

It is possible to classify perturbations using properties 
under the 4-dimensional coordinate transformation. Let 
us discuss the scalar, vector, and tensor perturbations, 
separately. 

A. Scalar Perturbation 

Although we take the Gaussian normal coordinate sys- 
tem for discussing the effective theory on the brane, the 
following gauge is convenient for the calculation in the 
case of the scalar perturbations: 



ds 2 



(1 + 6ip + 2^)dr 2 + (1 + 5<p + 2<S>)rj^dx fl dx v 



+r 2 (l-36<p)d9 2 



(21) 
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where we absorbed the background value of the deficit 
angle into the definition of 8, i.e. < 6 < 2ttc. In 
linearized Einstein Gauss-Bonnet gravity, we do not re- 
quire the GB term in the flat bulk background because 
it consists of a quadratic form of curvature. Then, the 
6-dimensional Einstein equation becomes 



-D5u> + -Sip" - —Sip' + - + 4$" 
2 2 2r r 



3$' + Sip 



1 

_ i 

r 
3 







OStp + -Sip" + —Sip' 
2r 



2 r 2 
f + 2$ 



1 , 4 , 

r r 



= , 



,(22) 
(23) 

(24) 
(25) 



1 



1 



-DSip + -Sip" + —Sip' + □$ + + = . (26) 
2 2 2r r 

Integrating Eqs. (|23|l and (|25|l . we can write 'J and <5<^ in 
terms of 



* = -2$ , 

3 r 
<5<£ = -r<f>' + $ + 2^( r ) ' 



(27) 
(28) 



where d(r) is a constant of integration and depends only 
on r. Here, we imposed the regularity at x^ — > oo. Af- 
ter eliminating using Eqs. |H|) and putting 
Eqs. (|27|l and (|28l) in the resulting equation, we find the 
equation of motion for 



-$' + □$ = 

r 



(29) 



After applying the same procedure to Eq. I|22|l and (|26|l . 
we obtain 



$" + _$' + □$+ [d( r )r]' = . (30) 
r 3r 

Comparing the above Eq. (|30() with Eq. (|29|l , we get 



d(r)r — const. = Co 



(31) 



This constant is related to the mass of the system deter- 
mined from the given energy-momentum tensor Q . 
If we define Fourier transformation, 



(32) 



*(r,a^) = / d 4 pe ip - x <£{r,p) 
Then the solution of Eq. l(2"9"j) is 

$(r,p) - — Ji(gr) , 



(33) 



where = — and we have imposed the regularity 
at the origin. The amplitude A(q) should be determined 
by the matching condition. The solution is given by 

ds 2 = (1 + 8g rr ) dr 2 + (?/ M „ + Sg^) dx^dx" 

d9 2 , (34) 



2 r 2 



(35) 



where 
Sg rr 

Sg^u 

Sgee 

Here, we have defined 

Zx(t,x») = J d 4 pe ip x A(q)J 1 (qr) , (38) 
Z2(r,aT) = ( d i pe tpx A(q)qJ 2 (qr) . (39) 



_ Co --Z 1 (r,x») + -Z 2 (r,x> 1 )) (37) 



In order to impose the matching conditions, we must 
move to the Gaussian normal coordinate system 



ds 2 



dr 2 + (r)^ + %,) dx»dx v + (1 + Sgee) d9 2 , (40) 



where the axial symmetry is taken into account. This 
can be achieved by the gauge transformations associated 
with the infinitesimal coordinate transformations x^ — > 
x* -f: 



where 



<5g ee = <5gee + 2c 2 r£ r 



dr5g r 



ZfJ, = o / dr / drS 9rr,n ■ 

1 Jo Jo 



(41) 
(42) 

(43) 
(44) 



Note that the location of the brane is not necessarily at 
constant r. However, there exist residual gauge transfor- 
mations 



Zr 



(45) 
(46) 



Using this residual gauge, one can adjust the brane posi- 
tion to be located at the constant r in the new coordinate 
system. Thus, the distortion function x( xfJ ") takes into 
account the fact that the brane is not necessarily located 
at the constant r in the coordinate system which we used 
to solve equations in the bulk. The final result is written 
by 



ds 2 = dr 2 + {\ + h^)i 1 ^dx^dx' J +L 2 d9 2 , (47) 



where 
hiijj 



C± + l Zl {r,x»)-lz 2 {r,x») 



dr 



dv8Qj- r in/ 

V 1 1 



(48) 



2 2 

c r 



Co--Z 1 {r : x») + lz 2 {r,x») 
r I 



-c r drSg rr + 2c r\ 
lo 



(49) 
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and we have reintroduced the background deficit, c, ex- 
plicitly. Now, the brane is located at r = in this coor- 
dinate system. 



B. Vector perturbations 

Now, we move on to vector perturbations which turn 
out to be irrelevant. As to the vector perturbations, we 
choose the gauge 

ds 2 = dr 2 + { Vltl/ + F^ v + F Vjfl ) dx»dx v + c 2 r 2 d9 2 ,(50) 

where F ^ ^ = 0. Note that there exist residual gauge 
transformations 

F M = F M + ^«) , = 0, (51) 

where the parameter depends only on x^ . 
The equation of motion for F^ becomes 

(F^u + Fv,n)" + \ [F^ v + F„J = . (52) 

The regular solution of this equation is 

F„ = /„(*") . (53) 

This solution can be eliminated using the residual gauge 
transformations Q51JI. Hence, there exist no vector per- 
turbations. 



C. Tensor perturbations 

Finally, we consider tensor perturbations. They are 
characterized by a metric 

ds 2 = dr 2 + (n^ + dx»dx v + c 2 r 2 d6 2 , (54) 

where hJJ satisfies transverse and traceless conditions. 
Tensor perturbations remain invariant under a coordi- 
nate transformation. The equation of motion for /iJJ 
gives 



<r + ^ T + n^' = 0. (55) 
The solution becomes 

h™{r,x») = J d A p e^e, u (q)Mqr) , (56) 

where should be determined by the matching condi- 
tion. 

V. EFFECTIVE THEORY 

A. Linear Regime 

Since we have solved the bulk geometry, we next 
discuss the matching conditions. Using the energy- 
momentum tensor (|12[) . matching conditions Eqs. (JSJ) and 



© reduce to 

-8K' = n 2 S 9 e S(r) , 

SK'i* v - 6%5K' = K 2 S^J(r) 



(57) 
(58) 



SH6L" - AaV'SG^u = — T"„<S(r) , (59) 

where we dropped all the remaining terms in Eqs. l|8"jl 
and @, because they do not contribute at first order in 
perturbation. 

Using Eq. (|49fl . the deficit angle can be obtained from 



lim L'(r = e) = Um \b'(r = e) + 5L'(r = e)l 
= c(l - C ) . 



(60) 



Note that the deficit angle is perturbed on the brane at 
the linear level, but it is constant. From Eqs. (|48|) and 
(|56|l . the extrinsic curvature in this limit gives 



lim Kuu(r = e,x^) 



lim 5 K, 



Kmti (r = e) = x,nu (61) 



Notice that the background value of is zero and the 
tensor perturbations do not contribute to the extrinsic 
curvature (|61|l . 

From this result, we see that the extrinsic curva- 
ture is nonzero near the brane, if the distortion field 
x(x M ) exists. In the coordinate system where we have 
solved the equations of motion in the bulk, the loca- 
tion of the brane is specified by x( xfl ) an d the result- 
ing shape of the braneworld becomes the deformed cylin- 
der in which the distortion is represented by x( x ^) ( see 
Fig.l). At first sight, this deformed cylinder looks like 
a 5-dimensional hypersurface. However, at the location 
of the brane (r = e in the Gaussian normal gauge), the 
circumference radius of the cylinder vanishes because of 
limL(r = e, x^) = there. Hence, the deformed cylinder 

is the 4-dimensional braneworld. 

Now, we must specify K^ u (r = 0,^) which is not 
known a priori. Without losing the generality, we can 
write 

K„ v {r = 0, x») = xoM^) + M^) . P% = >( 62 ) 

where we decomposed the extrinsic curvature at r = 
into the traceless part p^{x^) and the trace part 
Dxoix^). As the field x( xfi ) represents the location of 
the brane r = e, it is natural to regard Xoi^) as the 
location of r = 0. Hence, ^(x^) = x( xfi ) — Xo(^ M ) can 
be interpreted as the thickness of the braneworld. 

With these considerations, we find Eqs. JS7J, l|58|) and 
(|5*97l become 

n( x -Xo) = -K 2 S e e , (63) 

r% - %n x - [xo'% - ssnxo] - p\ = * 2 s\ ,(64) 

Si [8ti - SL' Q ] - 4a [b' - b' ] 5G» V = ^-T\ . (65) 

Z7T 
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L(r=e) =0 




FIG. 1: Schematic picture for the quasi-thick braneworld. 
The deformed cylinder has zero circumference radius. Hence, 
it is a 4-dimensional spacetime. ip^x^) = x{ x>1 ) ~ Xo{ x>1 ) 
represents the thickness of the brane. 



Adopting the standard boundary condition L'(r = 0) = 
b'(r = 0) + 5L'(r = 0) = 1 even in the perturbed space- 
time, Eq. H65JI leads to 



Lx v 



87ra(l - c) 



cC 



*4a(l-c) ' 



(66) 



where the second term of the right hand side of (|66|) is 
the cosmological constant. The Einstein tensor includes 
the contribution from the tensor perturbations. Eq. I|66|) 
proves that Einstein gravity is recovered for the codimen- 
sion 2 braneworld in the case of the linearized gravity. 
This fact is consistent with the result in 0- However, 
this is not the end of our story. We must also solve the 
conditions l|63l) and (|64|l . 

The conservation law T a e - a = for the e = /x at the 
linearized level gives the 4-dimensional conservation law, 











(67) 



where | denotes the 4-dimensional covariant derivative. 
On the other hand, the conservation law for the e = r 
implies, 



SK ^- + b'S 6 



6(r) = 



(68) 



Note that Eq. (|68|l appears to contain an ambiguity, since 
K^ v must be evaluated at r = where it is discontinuous. 
Thus we define the following quantities 



8K = lim - [SK(r = e) + SK(r = 0)] 

e^0 2 



]im-[b'(r = e) + b'(r = 0)] 

£— >0 2 



Using these quantities, one can deduce 
a 



SK^- + b'S e 







(69) 
(70) 

(71) 



The conservation law (|71|l gives the constraint on the 
distributional sources as 



1 



□(X + Xo) = -y 



-K 2 S 6 



(72) 



where we used Eq. Combining Eq. JB3J| with Uty . 

we find 



□x 



1 



-s° 



□Xo 



1 



(73) 



The above Eqs. (|73|l tell us that \a — C X- Since c < 1, 
ip = X ~ Xo = (1 — c)x is positive. Hence, it is legitimate 
to interpret ip as the effective thickness. In the meantime, 
comparing Eq. I|63|) with the trace of Eq. (|64|) , we find the 
relation between the components of energy-momentum 
tensor as 



s% = is" 



(74) 



Hence, once we give the 4-dimensional matter S^ Vl the 
extra component See can be determined. It should be 
noted that the effective thickness ip satisfies 



-5", 



(75) 



This is reminiscent of the radion in the codimension 1 
braneworld |lf)j | . The remaining quantity /5 M „ can be also 
obtained from Eq. (|64l) . 



B. Nonlinear Regime 

Our main point is the formulation of the framework of 
the quasi-thick braneworld. The importance of this point 
can be manifest in the next order calculations. Up to the 
second order, we expect 



87ta(l 
+c 



(2) 

V 



X'%X' Q 



cC ) 4a(l 
v - DXX'% 



■cCo) 



-6tt 



1 - c 



(76) 



where the corrections in the last three lines arises due to 

(2) 

W ^ and carry the information of thickness of the brane. 

(2) 

Another correction come from L ' carries the effects of 
bulk, that is, the deformed deficit angle. This needs the 
next order analysis. Thus, it turns out that some correc- 
tions due to the quasi-thickness can be expected at the 
second order. 

Here, we should point out a possibility that the above 
corrections to Einstein equations could mimic the ob- 
served features of dark energy/dark matter. To give some 



7 



flavor, let us extrapolate the effective theory l|76(l to the 
nonlinear regime. In the case of cosmology, one can put 



-PD 

pv5) 



(77) 



Then, one finds when the condition 



5k 4 



3(1 - c) 



Pd < 



(2) 

L 



cC 



4a(l 



cC ) 



(78) 



is satisfied, these kinds of corrections behave as dark en- 
ergy. While, when the condition 



(2) 

V 



cCn 



llK 4 



4a(l-c + cC* ) 9(1 -c) 



pI 



(79) 



is satisfied, the effective energy-momentum tensor due to 
corrections can be regarded as dark matter on the brane. 
To verify these expectations, we need to solve the bulk 
geometry. 

In reality, we need to consider the full nonlinear theory 
to discuss the cosmology. The nonlinear effective equa- 
tion should read 



lim 



SSL' - \u\l'G» v + L'W V \ 



4a 



6iU KUL-K a & L^ v 



)} 



— ,(80) 



where the left hand side Eq. (|8U|) represents the discon- 
tinuity at r = 0. To fully solve the system consistently, 
we must also take into account the matching conditions 
come from the less singular part which determines the 
boundary conditions completely. 

It is also intriguing to see how the gravitational waves 
on the brane are affected by the bulk geometry in the 
nonlinear situations such as the radiation from the parti- 
cle falling into the black hole. For this purpose, we must 
solve the above non-linear problem. 



VI. CONCLUSION 

We have studied the codimension 2 Einstein Gauss- 
Bonnet brancworld in the axially symmetric 6- 



dimensional spacetime. We have carefully examined the 
structure of possible singularities in the equations of mo- 
tion which characterize the braneworld through matching 
conditions. It turned out that the thickness of the brane 
can be taken into account within the context of the dis- 
tributional source, which we dubbed quasi-thickness. In 
the case of the linearized gravity, we have solved all of the 
equations of motion in the bulk and shown the conven- 
tional Einstein gravity can be recovered on the brane. In 
this process, all of the necessary boundary conditions are 
clarified. In the nonlinear regime, we found corrections 
due to the thickness and the bulk geometry. We stressed 
that the interplay between the bulk and the brane has 
been completely determined in the context of the quasi- 
thick braneworld. 

In our formulation of the codimension 2 braneworld, 
T^ v can be identified with the ordinary matter. While 
Sfj, v is a kind of dark energy/dark matter. The possibility 
that the thickness plays a role of the dark sector in the 
universe is attractive. It would be also interesting to 
construct a viable particle physics model in this context. 

There are many remaining issues to be solved. Once 
Tpv and S^ v are given, we can analyze the nonlinear grav- 
ity to reveal the nature of the corrections. This would 
be very complicated but the analysis is straightforward. 
Application to cosmology seems to be trivial in the linear 
regime. This is in contrast to the case of codimension 1 
braneworld where the analysis of the cosmological per- 
turbations are very complicated 11]. It is interesting to 
consider the black holes in the context of the codimension 
2 braneworld. The role of the Gregory-Laflamme insta- 
bility should be clarified and the property of the gravita- 
tional waves from the black hole must be analyzed [l2| . 
It is also intriguing to extend our formalism to the higher 
codimension braneworld Hal. 
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